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ABSTRACT 

While usually cosmological initial conditions are assumed to be Gaussian, inflationary theo- 
ries can predict a certain amount of primordial non-Gaussianity which can have an impact on 
the statistical properties of the lensing observables. In order to evaluate this effect, we build 
a large set of realistic maps of different lensing quantities starting from light-cones extracted 
from large dark-matter only N-body simulations with initial conditions corresponding to dif- 
ferent levels of primordial local non-Gaussianity strength /nl- Considering various statistical 
quantities (PDF, power spectrum, shear in aperture, skewness and bispectrum) we find that 
the effect produced by the presence of primordial non-Gaussianity is relatively small, being 
of the order of few per cent for values of |/nl | compatible with the present CMB constraints 
and reaching at most 15-20 per cent for the most extreme cases with |/nl| = 1000. We also 
discuss the degeneracy of this effect with the uncertainties due to the power spectrum normal- 
ization (78, finding that an error in the determination of as of about 3 per cent gives differences 
comparable with non-Gaussian models having /nl == ±1000. These results suggest that the 
possible presence of an amount of primordial non-Gaussianity corresponding to |/nl| = 100 
is not hampering a robust determination of the main cosmological parameters in present and 
future weak lensing surveys, while a positive detection of deviations from the Gaussian hy- 
pothesis is possible only breaking the degeneracy with other cosmological parameters and 
using data from deep surveys covering a large fraction of the sky. 

Key words: cosmology: theory - gravitational lensing: weak - cosmological parameters - 
large-scale structure of the Universe - Methods: N-body simulations 



1 INTRODUCTION 

In recent years, the interest for an accurate measurement of the 
amount of non-Gaussianity present in the primordial density field 
has largely increased. The main reason is that this test is now con- 
sidered not only a general probe of the inflationary paradigm, but 
also a powerful tool to constrain the plethora of its different vari- 
ants. Only the most standard slow-rolling models based on a single 
field produce in fact almost uncorrected fluctuations, which is the 
motivation of the common assumption (and large simplification) 
that their distribution is Gaussian. In general, small deviations from 
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Gaussianity are predicted even for the simplest inflationary models, 
while non-standard models, like the scenarios based on the curva- 
ton, on the inhomogeneous reheating and on the Dira c-Bom-Infeld 
inflation allow much more significant departures (see lSartolo et al.l 
12004 . and references therein). 



It has become common to quantify the level of primordial 
non-Gaussianity adopting the dimensionless non-linearity parame- 
ter /nl (see, e.g.,|Salopek & Bondll990l : lGanguill994l : IVerde et all 
uOOd : iKomatsu & Spergel 2001), that measures the importance of 
the quadratic term in a sort of Taylor expansion of the gauge- 
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invariant Bardeen potential ^: 
$ = ^L + /n 



(*L 



m)) 



(1) 



here $l represents a Gaussian random field. Hereafter we will 
adopt the so-called large-scale structure (LSS) convention, where 
$ is linearly extrapolated to the present epoclQ Moreover we 
will consider only the so-called local shape for non-Gaussianity, 
in which the bispectrum signal is larger on squeezed triangle con- 
figurations. For mor e details about other possible shapes we refer to 
iBartolo et all ( l2004 ); l\^rde & Matarresel ( l2009h : lBartolo & Riottol 
( 12009 1) and references therein. 

At present the most stringent constraints on /nl come from 
the cosmic microwave background (CMB) data. Their discrimi- 
nating power derives from the fact that its temperature fluctua- 
tions trace the density perturbations before the gravitational non- 
linearities modify their original distribution. Whatever is the spe- 
cific test adopted, all CMB analysis consistently allow only very 
small deviations f rom Gaussianity: for instance, analyzing the re- 
cent WMAP data. lKomatsu et alj J201(t) fou nd that /nl varies be- 
tween -13 and 97. while ISmith et all bOOsI) found -5 < /nl < 
104 ( at 95 per cent confidence level); see also iKomatsu et al.l 
( l2009h . 

Alternative and complementary constraints on /nl can in 
prin ciple be der i ved analyzing the LSS (for recent reviews, see, 
e.g. IVerdell201CI ; iDesiacaues & SelialduOlOn . As already evident 
from the first gener a tion of non-Gaussian N - body simulations 
JMessina et alj Il990l ; iMoscardini et alj Il99ll : [Weinberg & Cold 
Il992h . the presence of a positive (negative) skewness in the PDF 
of the primordial density field tends to favor (disfavor) the forma- 
tion of cosmic structures, inducing a different timing in the whole 
process of gravitational instability. However, to be fully exploited, 
this approach needs to be complemented with reliable methods to 
disentangle from the primordial signal the non-Gaussian features 
introduced by the late non-linear evolution and by the possible pres- 
ence of a non-linear bias factor. While lensing statistics entirely 
avoid the latter problem, to attack the former one it is necessary 
to make use of both analytic techniques (like high-order perturba- 
tion theory and the Time-Renormalization Group approach; see, 
e.g.. lTaruva et a l."2008'; 'Bartolo et al."2009S and full N-body sim- 
ulations (see, e .g., Kang et al. 2007; Grossi et al. 2007; Dalai et al. 
2008l;IViel et al. 2009i; lDesiacques et alj2009l ;|Pillepich et alj2009l ; 



Grossi et al.,,2009.) to properly calibrate the theoretical predictions 



in the non-linear regime. The large amount of theoretical work 
recently done in this direction allowed to better understand what 
is the size of the effects on the abundance of non-linear struc- 
turesjJMatarrese et al. 2000; Verde et al. 2000; Mathis et al. 2004; 



Kang et al^ ,2007,; .Grossi et aU i2007i. ,2009. ; Maggiore & Riotto 
2009l;lRoncarelli et al.l2009l). on the halo biasing ( Dalai et al. 2008; 



McDonald 120081 ; iFedeli et al.' 2009; Carbone et al. 2010), on the 



galaxy bisp ectrum fSefusatti & Komatsu 2007 ; Jeong & Komatsu 
I2OO9I; iNis himichi et al. 2009), o n the mass density dis tribution 
JGrossi et al. 2008; Lam & ShettJ '20091; 'Lam et alj|2009h . on the 
topology (Matsubara 2003; Hikage et al. 2008 ), on the integrate d 
Sachs-Wolfe effect ( Afshordi & Tollev 2008; C arbone etalj2008h 
on the Ly-Q flux from low-density in tergalac t ic me d ium JViel et al 



l2009t) . on the 21 cm fluctuations JCooravl I2OO6I ; IPillepich et al. 



^ We recall that on scales smaller than the Hubble radius <E> con'esponds to 
the usual Newtonian peculiar potential (but with changed sign). 
^ With the cosmological parameters adopted in this paper, this con'esponds 
to values for /nl larger by a factor of ai 1.3 with respect to the so-called 
CMB convention, where $ is instead extrapolated at 2 = 00. 



I2007D and on the reionization process JCrociani et al.ll2009f) . The 
first attempts of an a pplication to real ob servational data gave very 
encouraging results: Isiosar et al.l J2008 1. combining the bias mea- 
surements for two sa mples of luminous red gala xies and quasars, 
found /nl = 48+74; I Afshordi & TollevI J2008h , studying the in- 
tegrated Sachs-Wolfe effect (ISW) in the NVSS survey, derived 
.fNL = 354 ± 165; all error bars are at 2-a levejj. Very recently 
IXia et alj J2010l) found 32 < /nl < 152 from the analysis of the 
auto-correlation of the brightest NVSS sources on angular scales of 
several degrees. 

In this paper we will focus on estimating the weak lensing sig- 
nals in scenarios with primordial non-Gaussianity. Being based on 
the measurement of the shear effect produced by the intervening 
large-scale structure of the universe on the images of background 
galaxies, gravitational lensing is a direct probe of the total mat- 
ter distribution. For this reason it is considered one of the most 
powerful tools to constrain the main cosmological parameters and 
many dedicated projects are in progress or under study. A very ex- 
citing perspective is certainly op ened by the ESA Cosmic Vision 
project EUCLID jLaureiisll2009l) . currently under study: the goal 
of its wide survey is to obtain the shear measurements for about 
40 galaxies per arcmin^ on the entire extragalactic sky with galac- 
tic latitude fe > 30 (approximately 20000 deg^). The possibility 
of using the weak lensing signals to constrain also the amount of 
primord ial non-Gaussianity has bee n already explored by different 
authors. lAmara & Refregiea ( 120040 used a generalized halo model 
to study the impact on the estimates of the power spectrum nor- 
malization (T8 of primordial non-Gaussianity, modeled assuming 
various lognormal dis tributions for the density field. More recently, 
iFedeli & Moscardinil (2010) computed the power spectrum of the 
weak cosmic shear for non-Gaussian models with different values 
of /nl- In particular, they improved the halo model including more 
accurate prescriptions for its ingredients (mass functions, bias and 
halo profile), calibrated on the last generation of non-Gaussian N- 
body simulations. The application of this model to a survey hav- 
ing the expected characteristics of the EUCLID project showed the 
possibility of a significant detection of non-Gaussianity at the level 
of I/nl I ~ few tens, once the remaining parameters are held fixed. 

In this paper we investigate weak lensing statistics in non- 
Gaussian scenarios using numerical rather than analytical tools. 
Specifically, we will create weak lensing maps performing ray- 
tracing simulations through very deep light-cones extracted from 
high-resolution N-body simulations. The advantage of this ap- 
proach is twofold. First of all, N-body simulations permit to fully 
account for the non-linear evolution which is usually modeled less 
accurately by analytical means. Second, numerical experiments al- 
low us to extract a large set of realistic weak lensing maps that can 
be used for better evaluating the statistical robustness of the results. 

The main goal of our numerical work is to figure out what 
are the observational evidences of the presence of some level of 
primordial non-Gaussianity, as quantified by the /nl parameter. In 
particular we will compute a large set of weak lensing statistics in 
models with different /nl and we will quantify the deviations from 
the corresponding results in the Gaussian scenario. This is impor- 
tant not only to address the possibility of a positive detection with 
future data, but also to establish at which level an amount of pri- 
mordial non-Gaussianity compatible with the present observational 



^ We report the values as revised bv lGrossi et alj j200 y) to include a cor- 
rection mimicking the ellipsoidal collapse and converted to the LSS con- 
vention. 
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constraints can hamper an accurate measurement of the other cos- 
mological parameters. We must recall that Gaussian initial condi- 
tions are virtually always assumed in their practical derivation. 

The plan of this paper is as follows. In Section |2] we review 
the basis of the lensing formalism necessary to the present work. 
In Section [3] we describe the cosmological simulations and the 
numerical procedure to build the lensing maps. In Section |4] we 
present our main results about the statistical properties of the dif- 
ferent lensing quantities investigated: the probability distribution 
function, the third-order moment, the power spectrum and the bis- 
pectrum. In Section[5]we compare the effects produced by primor- 
dial non-Gaussianity to the uncertainties related to power spectrum 
normalization as- Finally, in Sect. |6] we draw our conclusions. 



2 LENSING THEORY 

In this section we give a short summary of the aspects of the the- 
ory of gravitational lensing that will be u sed throughout this work. 
For m ore detail we refer to the review bv lSartelmann & Schneiden 
1 I2OOII) . 

We can start by describing how light rays are deflected by 
the presence of structures in the universe. Since the coherence 
length of cosmic structures is small compared to the Hubble ra- 
dius, it is possible to slice the large-scale structures into shells and 
use the so-called thin-screen approximation for them, which al- 
lows to consider only the projected mass distribution of the slices. 
Denoting by E(6') the projected mass distribution of the lens at 
the angular position vector 6, the convergence can be defined as 
K{e) = E(^)/Ecrit, where 



Dds 



A-kG DaDs 



(2) 



is the critical surface density. In the previous equation Ds, Dd and 
Dds represent the angular-diameter distances between the observer 
and the source, between the observer and the lens, and between the 
lens and the source, respectively. 

Thanks to the thin-screen approximation, the object acting as 
lens can be completely described using its lensing potential '^, 
which is related to the convergence through the two-dimensional 
Poisson equation 



V^*(6') = 2k{9) 



(3) 



The deflection angle a is simply the gradient of the lensing poten- 
tial, a = V*. 

Up to second order, lensing-induced image distortions are 
given by 



AijOj + —DijkOjOk 



(4) 



(see, e.g., iGo l dberg & B acorJ l2005l : iBacon et al.l |2006|) . where 
Aij = djO'i represent the elements of the Jacobian matrix of the 
lens equation, 6/ are the unlensed coordinates, and the tensor is 
defined as Dijk = dkAij. In the previous equations di = d/d9i. 
The quantities A and D can be conveniently expressed as a func- 
tion of the convergence «:, of the complex shear term 7 = 71 + J72, 



and of their derivatives: 

1 - K - 71 

-72 
-271,1 —72,2 



A 



D 



ijl 



-72,1 



-72 
1 - K + 71 

-72,1 
"72,2 



(5) 



D 



ij2 



-72,1 
-72,2 



—72,2 
271,2 — 72,1 



in the previous equations the comma indicates the derivative. The 
shear derivatives can be combined to construct two new quantities, 
F and G, called ^'rit and second flexion respectively, defined as 

F = Fi+iF2 = (71,1 + 72,2) + 4(72,1 - 71,2) (6) 

G = Gi +iG2 = (71,1 - 72,2) +2(72,1 +71,2)- (7) 

The previous formalism can be easily generalized to the case 
in which a continuous distribution of matter is considered. The vol- 
ume between the observer and the sources can be divided in a se- 
quence of sub-volumes having a size along the line-of-sight suffi- 
ciently small compared to the distances between the observer and 
the sub-volumes, and between those and the sources. The matter 
of each sub-volume can be projected onto a plane and then we 
are allowed to use again the thin-screen approximation described 
above. The final quantities, estimated on the source plane, will be 
the weighted sum of the relevant quantities, where the weight is 
given by a suitable ratio of the involved distances. In multiple lens- 
plane theory, rotation of light bundles can occur in addition to shear 
and convergence, but numerica l simulations have shown that these 
are negligibly small ( Jain et alJ uOOO). Then, all the lens properties 
are contained in the lensing potential. Knowing it on all lens planes 
allows us to write recursion relations representing the deflection an- 
gle, the shear, the effective convergence and the two flexions on the 
source plane. These relations will be given and discussed in more 
detail in Section[331 



3 THE LENSING SIMULATIONS 

3.1 The cosmological simulations 

To study the effect of non-Gaussianity on the weak lensing statis- 
tics, we use the outputs of a set of N-body cosmological simulations 
following the evolution of dark matter particles only. These simula- 
tions have been alre a dy use d bv lGrossi et alj feOOTn.lHikage et al.l 
(2008). Gr ossi et al.l J200^ and (Roncarelli etalj ( 12009 V Here we 
will summarize the information relevant for this paper, referring 
the interested readers to the original papers for further details. 

The simulations were carried out using the publicly available 
code GADGET-2 ( Springe! 2005) assuming a ACDM model with 
total matter density ilm = 0.3, baryon density fib ~ 0.04 and cos- 
mological constant density SIa = 0.7; the Hubble parameter is set 
to h — 0.7, while the power spectrum is normalized to ag = 0.9. 
The simulated box has a comoving length of 500 Mpc/h and con- 
tains 800'^ particles with a mass resolution of 2 x 10^'^ Mq/H. 
The comoving softening scale is epi — 12.5 kpc/h. The different 
outputs are equispaced in comoving space by 250 Mpc/h. The set 
is composed by seven cosmological simulations considering differ- 
ent values of the non-Gaussianity parameter /nl: ±1000, ±500, 
±100, plus the standard Gaussian case /nl = 0. Notice that in this 
work we prefer to make use of this set of N-body sim ulations in 
spite of the more recent one presented in lOrossi et al.l J200ft) . The 
reason for this choice is twofold. First, the mass resolution is bet- 
ter and this allows us to have more robust results on small scales. 
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in particular for shear and flexions; second, the assumed values for 
the parameter /nl span a larger range, covering also cases where 
the weak lensing signal produced by primordial non-Gaussianity is 
larger. 

3.2 Building tlie mock light-cones 

To construct the mock light-cones used to perform the ray-tracing 
simulations, we follow the same procedure described in lPace et alJ 
( l2007r) . to which we refer for more details. Since the different out- 
puts of an N-body simulation represent the redshift evolution of the 
same initial matter distribution, we need to apply a specific pro- 
cedure to avoid the introduction of biases related to the fact that 
the same structures appear approximately at the same positions in 
different outputs. For this reason we randomly shift and rotate the 
particle positions exploiting the periodicity of boundaries. As said 
before, given our choice for the output redshifts, they overlap for 50 
per cent of their comoving side-length, therefore we can just con- 
sider particles in the lower half of the rotated and translated boxes. 

Once selected the particles to be used, we project them per- 
pendicularly to the line-of-sight on a regular two-dimensional grid 
and compute the projected mass density fiel d using the triangular- 
shape d-cloud (TSC) mass assignment dHocknev & Eastwood 
Il988l) . Finally via FFT techniques it is possible to recover the lens- 
ing potential associated with the considered matter distribution. 
More details on the numerical procedure will be given in the next 
section. 

Our mock light-cones extend along the line-of-sight up to red- 
shift Zs — 4, using A'out ~ 20 outputs. As all the distances are 
in comoving units, the opening angle of the ray-tracing simulation 
can be calculated using the comoving distance of the last plane in 
the stack, resulting in 6* = 5.876 degrees. Using a grid of 2048^ 
points, the corresponding angular resolution of the produced maps 
is 10.33 arcsec. 



3.3 Ray-tracing simulations 

Ray-tracing simulations consist in tracing back a bundle of light 
rays through the matter distribution of the light-cone, from the ob- 
server to the sources, which we place at redshift Zs = 4. 

Projected mass maps are converted to projected density- 
contrast maps 5{JJj°^'' by 






Ap 



Li 



(8) 



where Ai is the the area of the grid cell on the i-th plane, M}^ is 
the mass projected on the grid cell with indices {I, m) belonging to 
the i-th plane and Li = 250 Mpc/h is the depth of the i-th subvol- 
ume used to build mock cone. The density con trast is deliberately 
defined such as to have the unit of a length (see lHamana & Mellien 
uOOlh . In the previous equation p represents the average comoving 
density of the Universe. 

The lensing potential on each plane ^i is related to the pro- 
jected density contrast through the two-dimensional Poisson equa- 
tion, namely: 



V'*.(e"*) = Sfin, ( ^ 



'\0) ■ 



(9) 



Exploiting the periodic boundary conditions of the projected maps, 
Eq. ^ can be solved adopting the FFT techniques. Having the po- 
tential on each plane, the lensing quantities can be derived adopting 
standard finite-difference schemes. 



In order to perform the ray-tracing simulations, we need to ap- 
ply the multiple-plane theory to compute the total effect taking into 
account the contributions from each single lensing plane. A light- 
ray is deflected on each plane by the amount aiifii), thus the total 
deflection is given by the sum of all contributions. In particular, if 
the light-cone is sampled into A'out lens planes and the sources are 
located on the A^out + 1 plane, the relation giving the deflection 
angle on the i-th plane of a ray with image position 6\ reads (in 
comoving units) 



E 



fK{t^t 



■Wk 



}K{wi)ak 



Vs*fc(a 



(10) 



here ^^ (£) is the unsealed lensing potential, i.e. the Newtonian po- 
tential projected along the line-of-sight, /k is a function depending 
on the cosmology, w is the comoving distance and a represents the 
scale factor of the lens plane. Note that in general the light-rays 
will intercept the lens planes at arbitrary points, while the poten- 
tial is defined only at the grid points. Thus, it is necessary to use a 
bi-linear interpolation to compute the lensing quantities. 

Differentiating Eq. l llOb with respect to di and defining Ai = 
d6i/d6\ and Uu the matrix containing the second derivatives of the 
lensing potential, one obtains 



A. = I^Y. 



fK(Wk)fK{Wi - Wk) 

fK{wi)ak 



UkAk 



(11) 



In the previous equation, / represents the identity matrix. On the 
source plane, the lacobian matrix AjVo„t+i is given by 



A, 



-72 



-7i 

- OJ 



-72 +ljj 
- K + 71 



(12) 



where k is now the effective convergence and 7 = 71 + i72 is the 
effective shear. The term oj, called rotation, represents the asymme- 
try introduced by multiple lenses. 

Differentiating Eq. Jilt with respect to 9i, a recursive relation 
for the two flexions can be obtained: 



^l'^ = -E 



-T-j—T- IfkWkGu +UkD^' ] , 

(13) 
where Gu = Vjt/ is a tensor containing the third derivatives of 
the lensing potential. 

On the source plane, the tensor D reads as 



^JV+l 



D 



JV+l 



—271,1—72,2 -72,1+^1 

— 72,1 — Cdl ^72,2 

—72,1 —72,2 + OJ2 

—72,2 — OJ2 271,2 — 72,1 



(14) 



where, as for the Jacobian matrix, uji and 012 represent the asym- 
metric terms. In Eq. il4t we do not take into account the addi- 
tional terms, called twist and turn, introduced bv lBacon & Schafen 
(2009), but several tests performed at different resolutions assure 
us that their inclusion would not significantly affect our results. 



4 RESULTS 

In this section we present the results of our analysis of the statisti- 
cal properties of several lensing quantities extracted from the ray- 
tracing simulations described above. In particular we discuss the 
probability distribution function (PDF) in Sect. 14. 11 the third-order 
moment (skewness) in Sect. 14.21 the power spectrum in Sect. 14.31 
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the shear in aperture in Sect. 14.41 and finally the bispectrum in 
Sect. 1431 



4.1 Probability distribution functions 

We start considering the PDFs for several important lensing quan- 
tities. All our results are grouped in Fig. [T] The panels on the left 
present the actual PDFs of a given quantity, as obtained averag- 
ing over 60 different light-cone realizations. Given the smallness of 
the differences, we prefer to show only the results for the Gaussian 
model /nl = (red points and lines) and for the most extreme non- 
Gaussian models: /nl ~ —1000 (blue lines) and /nl ~ +1000 
(green lines): models with smaller non-Gaussianity are intermedi- 
ate between the displayed lines. For reference, we also show (for 
the Gaussian model only) the error bars, referring to the r.m.s. of 
the corresponding sample of realizations. Finally, the panels on the 
right present for all models (represented by different colour lines, 
as described in the label) the ratio between the non-Gaussian and 
Gaussian results. 

The top row of Fig. \T\ refers to the modulus of the deflec- 
tion angle a. The presence of primordial non-Gaussianity affects 
the PDFs by at most few percent in the most extreme cases with 
/nl ~ ±1000, while the models with primordial non-Gaussianity 
more consistent with current constraints (/nl = ±100) are almost 
indistinguishable from the Gaussian case. We also notice that the 
differences produced by non-Gaussianity are more evident in the 
high-value tails. These large deflections reflect lensing events asso- 
ciated with rare, large structures. In this sense these results confirm 
those on the halo abundance: unlikely events (large haloes) are use- 
ful probes of primordial non-Gaussianity. However, we have to note 
that in the high- value tails the statistics can be quite poor and noisy: 
a good modelization of the strong non-linear effects acting on the 
same scales is necessary to allow a positive detection. 

In the panels in the second row of Fig. [T] we show the analo- 
gous plots for the (modulus of the) shear 7. Since a negative value 
of the /nl parameter favors lower values of 7, we find that the ratio 
between non-Gaussian and Gaussian models is in this case larger 
than unity; the opposite trend holds for models with positive pri- 
mordial non-Gaussianity. As already pointed out for the deflection 
angle, the effect of a mild non-Gaussianity on the shear distribu- 
tion is tiny, below 1 per cent, which is comparable to the size of the 
error bars (shown only in the left panels). Only models with very 
high values of /nl display a deviation with respect to the Gaussian 
case that can be as large as 10 per cent. 

The PDFs for the first and second flexion are shown in the 
panels in the third and forth rows of Fig.[T] respectively. The trend 
is very similar to the case of the shear: compared to the Gaussian 
model, models with negative /nl have a higher probability of an 
excess at low values; the opposite trend holds in the high-value tail. 
The differences between Gaussian and non-Gaussian models are 
slightly more evident for the second flexion G. 

Finally the panels in the last row of Fig. [T] refer to the effec- 
tive convergence, for which we find a slightly different situation. In 
this case, in the maps there are also negative values, corresponding 
to underdense regions where the non-Gaussian PDFs show large 
deviations from the Gaussian results: for k — —0.2 we find 2-3 
per cent for /nl = ±100, « 15 per cent for /nl ~ ±500, and 
approximately 20 per cent for /nl = ±1000. No t e that this re- 
sult agrees with the analysis made by iGrossi et al.l ( 120081) on the 
three-dimensional density field and its probability distribution: un- 
derdense regions tend to maintain the imprint of primordial non- 
Gaussianity, suggesting that statistics based on voids can be a pow- 



erful tool to estimate /nl (see also lKamionkowski et al.l2009B . The 
plot shows that also the overdense regions keep the important im- 
prints of the primordial non-Gaussianity: for k — 0.5 we find 4-5 
per cent for /nl = ±100, ~ 25 per cent for /nl = ±500, and 
more than 60 per cent for /nl ~ ±1000. However, in this regime, 
the statistics is quite poor and the ratio between non-Gaussian and 
Gaussian PDFs becomes very noisy. As final comment on PDFs 
(and on their moments, see the following sections), it is important 
to stress that the lensing quantities here discussed are observables 
that directly measure the dark matter distribution. This is not true 
for the density PDF derived from galaxy surveys, where the com- 
plex effect of bias must be corrected for. 

4.2 Third-order moment 

A possible alternative way to detect the signatures produced by the 
presence of some primordial non-Gaussianity is to look at the high- 
order moments (skewness, kurtosis, etc.) of the distribution of the 
various lensing quantities. The power of this approach was already 
evident in the investigation of the corresponding qua ntities related 
to the density field (see, e.g.. lColes et alJll993l ; lGrossi et al.ll2008l; 
iLam & ShethluOOgh . However, the measurement of these statistics 
is often affected by large error bars, that increase with the order and 
hamper their reliable application to the real data. For this reason, 
in this section we will focus only on the skewness ^3 of the two- 
dimensional convergence field k^.j, defined as 



M3 = 



M2 



3/2 



iV2 



■I](«^'>. 



-\3 



(15) 



where k and fj,2 = X^i i{i^i,j~i^)^ /^^ represent the mean and the 
variance, computed on all A'^^ pixels of the maps. To derive these 
quantities, we first convolve the convergence map with a Gaussian 
filter of angular radius 9, then we subtract the mean and finally we 
normalize it appropriately. 

In the upper panel of Fig.|2]we present the results for the skew- 
ness of the effective convergence, as extracted from the Gaussian 
simulation. For reference we also show the results for the two most 
extreme non-Gaussian models. As expected, fis (which is com- 
puted as average over 60 different realizations) is a decreasing func- 
tion of the filtering radius 9: increasing the smoothing reduces the 
non-Gaussian features introduced by the non-linear evolution. Fur- 
thermore the error bars, representing the r.m.s. over the set of dif- 
ferent light-cone realizations, are slightly increasing with due to 
increasing shot noise, since increasing the smoothing radius we are 
averaging over a small er number of circ les. This behavior is analo- 
gous to what found bv ljain et alJOOOOl) in the determination of the 
skewness of the effective convergence. 

More interesting is the bottom panel of the same figure, where 
we show the ratio between the third-order moments computed in 
non-Gaussian and Gaussian simulations: we find a variation of 
the order of 15-20 per cent for /nl = ±1000, « 8 per cent 
for /nl = ±500 and only « 1 per cent for /nl ~ ±100. It 
is worth noticing that this effect is almost constant on scales up 
to 10 arcmin. Comparing these results with the size of the error 
bars shown in the upper panel, we can conclude that with relatively 
small maps, like those analyzed in this paper, we can use the skew- 
ness statistics to detect only strong primordial non-Gaussianities 
(/nl = ±500, ±1000). As shown by the results of the analysis 
made bv lFedeli & Moscardinil ( I201C1) . only with data covering very 
large area it would be possible to disentangle the effects of non- 
Gaussian models with /nl = ±100. Finally we notice that a sim- 
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Figure 1. Left panels: the differential PDF for different lensing quantities: the (modulus of the) deflection angle a (first row); the (modulus of the) shear 
7 (second row); the (moduli of the) two flexions F and G (third and fourth rows, respectively); the effective convergence k (last row). Points and error 
bars represent the average and rm.s. of a set of 60 different realizations for the Gaussian model, while blue and red lines present the average values of 
the simulation set corresponding to the two most extreme non-Gaussian models, /nl = —1000 and /nl = +1000 respectively. Right panels: the ratio 
between non-Gaussian and Gaussian differential PDFs for the same quantities shown in the left panels. Different color lines refer to models with different 
levels of primordial non-Gaussianity: /nl = +1000 (green), /nl = +500 (blue), /nl = +100 (cyan), /nl = —100 (orange), /nl = —500 (red) and 
/nl = -1000 (brown). 
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Figure 2. Upper panel: the skewness /13 of the effective convergence for 
the Gaussian model as a function of the angular' scale 6. The data and error 
bars represent the mean and the r.m.s. computed on a set of 60 different 
light-cone realizations, while blue and red lines present the average values 
of the simulation set corresponding to the two most extreme non-Gaussian 
models, /nl = —1000 and /nl = +1000 respectively. Lower panel: the 
ratio between the skewness results for non-Gaussian and Gaussian models. 
Different color lines refer to different values of /nl. as indicated in the 
labels. 
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Figure 3. Comparison between theoretical predictions and our numerical 
estimates, derived as average of 60 different realizations of the power spec- 
tra. Error bars correspond to the r.m.s. over the different realizations. Vari- 
ous lensing quantities are shown: effective convergence, shear and reduced 
shear' (upper panel); first and second flexions (bottom panel). The results 
refer to the Gaussian case. 



ilar analysis performed on the shear field, which is directly related 
to the convergence one, would provide comparable results. 



4.3 Power Spectrum 

An important theoretical quantity, directly related to the mat- 
ter power spectrum, is the effective convergence power spectrum 
Pnif), which is defined as the squared modulus of the Fourier 
transform of k, averaged on the modes having a given multipole 
I. Starting from Pk(^), it is possible to derive analytic expressions 
for the power spectra of the other lensing quantities, such as those 
for the shear (P-y) and the two flexions (Pf and Pc), namely: 

p-,{e.) = p^e) (16) 

PfW = PG(e)^fP4£). (17) 

In Fig. [3] we show, for the Gaussian model only, the compari- 
son between the power spectra extracted from our simulations and 
the corresponding theoretical predictions including the non-linear 
effect. We consider different lensing quantities: effective conver- 
gence, shear and reduced shear g = 7/(1 — k) (upper panel) and 
two components of the flexion (bottom panel). We notice that the 
agreement between the numerical results and the theoretical expec- 
tations holds, as expected, only in a given range of wavenumbers. 
On large scales (small £), the departure from theory is due to the 
low number of modes we can use to average the numerical power 
spectra. On small scales (large £), the disagreement is produced 



by the mass resolution (particularly relevant for the flexion power 
spectrum) and by numerical artifacts produced by the weighting 
function appearing in the definition of a given power spectrum. We 
can use Fig. [3] as a reference to define the range of i on which the 
numerical results can be trusted: typically 10^ ^ ^ ^ 10*. 

In Fig.|4]we show the ratio between the power spectra derived 
from the non-Gaussian and the Gaussian simulations. We present 
here only the results for the effective convergence since the other 
power spectra can be obtained by introducing the suitable depen- 
dence on £, and then the corresponding ratios are identical. As ex- 
pected, the effects of the primordial non-Gaussianity become more 
evident when the absolute value of the parameter /nl increases. At 
very large scales, where the evolution is still in the linear regime, 
the Gaussian model and all the non-Gaussian ones here considered 
display the same power spectrum, and the ratio is in practice unity. 
At smaller scales (£ > 100 — 200) the behaviors start to be differ- 
ent: a positive (negative) /nl parameter implies a larger (smaller) 
power spectrum with respect to the Gaussian case. It is also worth 
noticing that at very small scales {£ > 1000 — 2000) the ratio tends 
to an asymptotic value: this is due to the fact that in this range the 
power spectrum can be very well approximated by a power-law re- 
lation and all models turn out to have approximately the same slope. 
We have to recall, however, that some numerical artifacts can affect 
our results on these small scales. Looking at the plot, we see that, 
as expected, the differences between Gaussian and non-Gaussian 
models are tiny, for the power spectrum too: for the most extreme 
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Figure 4. The ratio between non-Gaussian and Gaussian estimates for the 
convergence power spectrum (averaged over 60 different realizations) as a 
function of the multipole I. Different color lines refer to various values of 
/iML as indicated in the labels. 



models (/nl = ±1000), the deviations are of the order of 10-12 
per cent, for the models with /nl = ±500 they are approximately 
5-6 per cent, while for the models with /nl = ±100 the deviations 
are smaller than 2 per cent. 



4.4 Shear In aperture 

A quantity which is directly related to the shear power spectrum 
but which can be more easily measured in real data is the so-called 
shear in aperture 7av • It represents the variance of the shear field 
within a circular aperture of a given radius 6: 



Numerical i — i — i 

Analytical 

ff,L=+1000 

fM,=-1000 




Figure 5. The shear in aperture. Top panel: comparison between the theoret- 
ical analytical relation and our estimate, derived as average of 60 numerical 
realizations; error bars represent the rm.s. The results refer to the Gaussian 
model. Blue and brown lines present the corresponding results for the two 
most extreme non-Gaussian models, /nl = +1000 and /nl = —1000 
respectively. Bottom panel: the ratio between the results for non-Gaussian 
and Gaussian models: different color lines refer to various values of /nl ^s 
indicated in the labels. 
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(18) 



where Ji is the first-order Bessel function of the first kind. 

The top panel of Fig. |5] shows the comparison between the 
theoretical expectation and the results derived from our Gaussian 
simulations (averaged over 60 different light-cone realizations): the 
good agreement is evident up to scales of S ~ 20 arcmin. For refer- 
ence we also plot the results for /nl = ±1000. The bottom panel 
of the same figure presents the ratio between the results for non- 
Gaussian and Gaussian models as a function of the aperture angle. 
Being a weighted integral of the shear power spectrum, |7av(S)P 
has a trend which is quite similar to the one displayed by P^,: the 
larger the absolute value of the /nl parameter is, the larger are the 
deviations from the Gaussian expectations. An increment (decre- 
ment) in the signal of the order of 10 per cent is expected for the 
model with /nl = ±1000 (/nl = -1000), while the differences 
reduce to less than 1 per cent for the models with /nl = ±100. As 
expected, we find that the deviations from the Gaussian expecta- 
tions are larger for small apertures, but rapidly wash out for angles 
larger than few arcmin, since the average of the shear involves more 
and more large scale structures. 



4.5 Bispectrum 

In several works (see e.g.lVerde et alj|l998l . l2000l : iTakada & Jainl 
|2004 ISefusattiet"ai]|2006l l2010h it has been shown that a large 
amount of cosmological information can be inferred from the study 



of high-order spectra (the so-called poly-spectra). Here we just 
limit our investigation to the simplest high-order spectrum, the bis- 
pectrum. We recall that for a Gaussian random field all available 
information is contained in the power spectrum, while all higher- 
order spectra vanish or are combinations of the power spectrum. As 
a consequence the measurement of the bispectrum is quite sensitive 
to the non-Gaussian features of the quantity under investigation. 

In general, the bispectrum depends on a triangular configura- 
tion and its values are very sensitive to the particular configuration 
adopted. Since the complete determination of the bispectrum is ex- 
tremely expensive from a computational point of view, we will re- 
strict our analysis to equilateral triangles only. Notice that, given 
our assumption of local shape for the primordial non-Gaussianity, 
the non-Gaussian signal on squeezed configurations is expected to 
be larger. 

In general the bispectrum of a given quantity / is defined as: 

B(/l, ^2, 012) = (27r)2fc(/i23)(/(/l)/(/2)/(/3)> ■ (19) 

In the previous equation, / represents the Fourier transform of the 
field /. Since the wavenumbers at which the bispectrum is com- 
puted (£i, £2 and £3) must form a closed triangle in Fourier space, 
they are related by the constraint i?i23 = ^1 ± £2 ± ^3 = 0. Fi- 
nally 9i2 is the angle between £1 and (.2 and fixes, together with the 
triangle condition, £3. 

Similar to Fig. [3] the top panels of Fig. [6] show, for the Gaus- 
sian model only, the comparison between our numerical estimates 
and the theoretical bispectrum computed for the effective conver- 
gence. The results from the simulations represent the average over 
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Figure 6. Top panel: comparison between theoretical predictions and our 
numerical estimates, derived as average of 60 different realizations of the 
bispectmm for the effective convergence. Error bars represent the r.m.s. of 
the 60 different realizations. The results refer to the Gaussian model only. 
Bottom panel: the ratio between non-Gaussian and Gaussian estimates for 
the bispectrum as a function of the multipole £. Different color lines refer 
to various values of /nl as indicated in the labels. 



60 different light-cone realizations with the corresponding r.m.s. 
plotted as error bars. The theoretical quantity are derived from the 
matter density bispectrum BpT{ki, fe, 612), that can be computed 
according to the perturbation theory. For instance, in the case of the 
effective convergence k, the bispectrum can be obtained by using 
the Limber's appr oximation in an ana logous way as for the power 
spectrum (see e.g. lCoorav & HuluOOlh . namely: 



B^h, £2,612) 



Jo 



Bp 



G^w) 



dw 

/x(w) 

/i £2 

/if(w)' /k(w) 



(20) 



where /x(u^) is a function depending on the cosmologi- 
cal parameters and Gjw) contains the distance weight 
JBartelmann & Schneideii I2OOII) . For the other lensing quan- 
tities, similar relations hold. In order to quantify the effects of 
primordial non-Gaussianity, the panel on the bottom shows the 
bispectrum computed for the different non-Gaussian models, 
normalized to the corresponding result for the Gaussian one. 

As for the case of the power spectrum, we see that the agree- 
ment between the theoretical expectations and the numerical pre- 
dictions holds only for a limited range of angular scales, but in 
this case the range for which the estimate of bispectrum is robust 
is much more limited. Focusing on scales where the numerical re- 
sults can be considered reliable, we find that the contribution of 
non-Gaussianity is not very relevant, reaching up to 15 per cent 
in the very non-linear regime, but only for the most extreme non- 
Gaussian models (/nl ~ ±1000). For models with /nl in the 



range allowed by the most recent CMB constraints, we do not ex- 
pect a difference with respect to the Gaussian case larger than few 
per cent. 



5 THE DEGENERACY BETWEEN PRIMORDIAL 

NON-GAUSSIANITY AND THE POWER SPECTRUM 
NORMALIZATION 

As shown in the previous sections, the effects due to primordial 
non-Gaussianity on the various weak lensing quantities are small, 
in particular when we consider models with values of /nl in the 
range allowed by present data. To better investigate the possibility 
of their future detection, here we will discuss their degeneracy with 
the effects produced by the variation of other relevant cosmological 
parameters. In particular we will focus on the degeneracy between 
the amount of primordial non-Gaussianity and the power spectrum 
normalization as, which is certainly the most important source of 
uncertainty. 

Since the different timing in the structure formation induced 
by /nl modifies the power spectrum and its amplitude, non- 
Gaussianity can be confused with a Gaussian model with a differ- 
ent as parameter. In particular a positive (negative) /nl has effects 
similar to a higher (lower) power spectrum normalization. In this 
section we show how much the uncertainty on the value of the as 
normalization affects the possibility of detecting non-Gaussianity. 
For the sake of simplicity, we restrict our analysis to the conver- 
gence power spectrum, the shear in aperture and the effective con- 
vergence bispectrum. 

In Fig.|7]we show a comparison between the contribution of 
non-Gaussianity and the changes produced by a variation of the 
spectrum normalization. In particular we consider as in the range 
between 0.87 and 0.93: the size of this variation approximately cor- 
responds to the typical 2a error derived from the most recent CMB 
analysis. It is clear from the figure that the error on the normal- 
ization has effects more important than the introduction of a weak 
primordial non-Gaussianity. The two most extreme models produce 
deviations from the Gaussian predictions that are of the same level 
of uncertainty Acrg = 0.03. On the contrary, the models with /nl 
in the range allowed by current constraints have effects well within 
the range allowed by the present uncertainties, since their size is 
similar to a variation of erg « 0.01. This implies that significant 
constraints on the amount of primordial non-Gaussianity can be ob- 
tained only reducing the uncertainties on the power spectrum nor- 
malization. 

We also notice that a lower (higher) normalization of the 
power spectrum does not change the overall shape of the quantities 
investigated and its effect is almost constant on all scales, while 
non-Gaussianity acts mostly at small scales (represented by large £ 
for the power spectrum and the bispectrum and by small 9 for the 
shear variance). This means that in principle it could be possible to 
distinguish the two effects observing the spectra at different scales. 
Thus, non-Gaussianity changes the shape of the curves much more 
than variations in as do. 



6 CONCLUSIONS 

In this work, we used the outputs of N-body simulations to cre- 
ate a large set of realistic mock maps for several lensing quantities 
(deflection angle, effective convergence, shear and the two compo- 
nents of the flexion) in the framework of cosmological models with 
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Figure 7. Comparing the effects of different levels of primordial non- 
Gaussianity to the ones produced by different power spectrum normaliza- 
tions. The panels show various lensing quantities: effective convergence 
power spectrum (upper panel); shear variance (central panel); effective 
convergence bispectrum (lower panel). The results are presented as ratio 
with respect to the same quantity computed for the Gaussian model with 
fT8 = 0.9 (solid horizontal line) and refer to Gaussian models with different 
(Tg normalizations: 0.93 (dot-dashed line), 0.92 (dot with long dash Hne), 
0.91 (dotted line), 0.89 (double dots hne), 0.88 (dashed line with three dots), 
0.87 (short dashed line). For reference we also show the effects produced by 
different levels of primordial non-Gaussianity: /nl = +1000 (cyan line); 
/nl = +500 (blue line); /nl = +100 (red line); /nl = -100 (magenta 
line); /nl = —500 (orange line); /nl = —1000 (brown line). 



different amount of primordial non-Gaussianity, quantified using 
the dimensionless parameter /nl- In particular we considered sev- 
eral statistical properties (PDF, power spectrum, shear in aperture, 
skewness and bispectrum) and compared the results with the corre- 
sponding ones obtained in the Gaussian case. Our main results can 
be summarized as follows. 

• For all quantities here considered the eifect produced by the 
presence of primordial non-Gaussianity is relatively small, amount- 
ing to differences of few per cent for |/nl| = 100, 5-10 per cent 
for |/nl| = 500, and 15-20 per cent for |/nl| = 1000. These re- 
su lts are in good agreement w ith the analytic predictions presented 
in lFedeU & Moscardinil ( l2010h . 



• The largest effects are visible on small scales (i.e. for large 
multipoles I > 1000), where, however, also non-linearity can pro- 
duce strong effects which have to be accurately modeled. 

• The most promising statistical tests to search for imprints of 
primordial non-Gaussianity are the (convergence and shear) power 
spectra and the (convergence) bispectrum, thanks to the smaller size 
of their error bars at the relevant scales. 

• The differences of the various PDFs in both rare-event tails 
are also important, but their discriminating power is reduced by the 
poor statistics and by the high-level of noise. 

• We compared the effects produced by the primordial non- 
Gaussianity with the uncertainties due to the power spectrum nor- 
malization: an error in the determination of erg of about 3 per cent 
gives differences comparable with the non-Gaussian models with 
/nl = ±1000, while for more realistic non-Gaussian models with 
/nl = ±100 the effects is even larger than the one induced by 
A(7g = 0.01. 

As said, a significant covariance exists between primordial 
non-Gaussianity and fundamental cosmological parameters, espe- 
cially with (Tg. The smallness of non-Gaussian effects found in 
our analysis means that one can obtain a precise estimate of these 
parameters despite their covariance with /nl- As an example, ig- 
noring the possible presence of a primordial non-Gaussianity with 
/nl ~ 100, consistent with current observational constraints, 
would induce a mere 0.2 per cent uncertainty in the estimate of 
(Tg, i.e. well below the current l-cr error. The same argument, in 
reverse, is telling us that the search for non-Gaussian features in 
the weak-lensing statistics is a very challenging task. According to 
our results, it can be successfully completed only if errorbars can 
be significantly reduced. For the purpose of discriminating among 
competing model s the reduction needs n ot to b e dramatic. Indeed, 
as pointed out bv iFedeli & Moscardinil ilOV j). the fact that devi- 
ations from Gaussianity are small but systematics allows one to 
estimate /nl by adding the statistical information from different 
angular bins. 

A more significant reduction in the errorbars is required to 
break the degeneracy between non-Gaussianity with cosmological 
parameters, say erg. Indeed, as shown in Fig. [7] removing such de- 
generacy requires to compare information fr om a limited number 
of bin s at very different scales. As shown bv lFedeli & Moscardinil 
yOlOt), both tasks could be achieved with next-generation, aU-sky 
surveys. A good example is certainly represented by the ESA mis- 
sio n EUCLID (Laureiis 2009) . According to the analysis made 
bv iFedeli & Moscardinil j2010l) and based on analytic predictions 
for the power spectrum quantitatively confirmed in this numerical 
work, we expect that the quality and the quantity of the EUCLID 
data will allow to constrain /nl at the level of few tens, opening 
the possibility of discriminating between the various inflationary 
models. 
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